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An  experiment  has  been  conducted  to  determine  die  contribution  of  atmospheric  scatter  to  the  severity  of  the 
dazzle  experienced  by  a  human  under  illumination  from  a  visible  laser.  A  15  W  532  nm  laser  was  propa^ted  over 
a  380  m  outdoor  range  in  San  Antonio,  Texas,  over  nine  data  collection  sessions  spanning  June  and  July  2014. 
A  narrow  acceptance  angle  detector  was  used  to  measure  scattered  laser  radiation  within  the  laser  beam  at  different 
angles  from  its  axis.  Atmospheric  conditions  were  logged  via  a  local  weather  station,  and  air  quality  data  were 
taken  from  a  nearby  continuous  air  monitoring  station.  The  measured  laser  irradiance  data  showed  very  little 
variation  across  the  sessions  and  a  single  fitting  equation  was  derived  for  the  atmospheric  scatter  function.  With 
very  conservative  estimates  of  the  scatter  from  the  human  eye,  atmospheric  scatter  was  found  to  contribute  no 
more  than  5%  to  the  overall  veiling  luminance  across  the  scene  for  a  human  observer  experiencing  laser  eye  dazzle. 
It  was  concluded  that  atmospheric  scatter  does  not  make  a  significant  contribution  to  laser  eye  dazzle  for  short- 
range  laser  engagements  in  atmospheres  of  good  to  moderate  air  quality,  which  account  for  99-5%  of  conditions  in 
San  Antonio,  Texas. 

OCtS  codes;  (010.1310)  Atmospheric  scattering;  (140.7300)  Visible  lasers;  (140.3360)  Laser  safety  and  eye  protection;  (230.5820) 
Scattering  measurements;  (330.4060)  Vision  modeling;  (330.4595)  Optical  effects  on  vision. 
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1.  INTRODUCTION 

Laser  eye  dazzle  is  the  temporary  visual  obscuration  caused  by 
visible  wavelength  laser  beams  being  incident  upon  the  eye. 
Commercial  airline  pilots  are  regularly  subjected  to  such  dazzle 
through  malicious  laser  use  by  individuals  on  the  ground  [1], 
while  “laser  dazzlers"1  are  often  used  in  military  and  security 
applications  as  nonlethal  weapons  with  the  aim  of  warning 
or  determining  intent  [2]. 

Spatial  variations  in  the  refractive  index  of  the  ocular  media 
serve  to  scatter  light  from  intense  light  sources  and  spread  their 
image  across  the  visual  field  [3],  The  subjective  angular 
subtense  of  a  laser  source  is  determined  by  several  factors, 
but  primarily  by  the  laser  irradiance  incident  upon  the  eve 
(higher  irradiance  =  larger  dazzle  field)  and  the  ambient  lumi 
nance  (lower  light  level  =  larger  dazzle  field).  Other  contribut 
ing  factors  include  the  wavelength  of  the  laser  beam,  and  the 
age  and  pigmentation  of  the  observers  eye  [4]. 

In  everyday  outdoor  conditions  the  above  described 
“intraocular  scatter1  is  not  the  only  contributor  to  the  severity 


of  laser  eye  dazzle,  as  there  are  several  extraocular  (external  to 
the  eye)  scattering  components  with  the  ability  to  influence  the 
dazzle  experience.  Optical  elements  between  the  eye  and  the 
laser  source,  such  as  spectacles  [5]  or  a  windshield  [6],  are  ex 
amples  of  eKtraocuIar  scatter  components  which  can  affect  the 
severity  of  laser  eye  dazzle. 

This  study  concerns  the  contribution  of  another  extraocular 
scatter  source  the  atmosphere  through  which  the  laser  beam 
propagates.  An  intense  laser  beam  can  easily  be  viewed  off  axis 
(Fig.  1)  because  of  atmospheric  scattering.  In  this  image  light  is 
scattered  away  from  the  axis  of  the  laser  and  toward  the  viewer. 
When  viewing  the  laser  on  axis,  single  or  multiple  atmospheric 
scattering  processes  can  cause  light  to  be  directed  toward  the 
eye  from  angles  outside  of  the  laser  axis,  as  shown  schematically 
in  Fig.  2.  This  light  can  add  a  “corona"1  of  light  around  the 
central  bright  point  source  a  spatial  spreading  of  the  laser 
beam  due  to  the  eye  receiving  light  from  a  wader  angular 
subtense  than  it  would  have  done  in  the  absence  of  any  atmos 
pheric  scattering. 
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Fig.  2.  Schematic  of  atmospheric  scatter  showing  single  and  multi 
pie  scattering  processes  causing  light  to  be  directed  toward  the  eye 
from  angles  outside  of  the  laser  axis. 


Atmospheric  scatter  of  visible  wavelengths  is  primarily 
caused  by  fine  particles  (<2.5  |im  diameter)  suspended  in 
the  air  which  can  diffract,  refract,  or  phase  shift  light  [7], 
and  occurs  predominantly  in  the  forward  direction  (i.e.,  the 
direction  in  which  the  light  is  traveling).  Such  fine  aerosols, 
typically  originating  from  combustion  processes,  are  more 
efficient  at  scattering  and  more  plentiful  in  the  atmosphere  than 
coarse  particles  such  as  road  dust  and  therefore  dominate 
scattering  processes.  Air  molecules  are  also  capable  of  scattering 
but  they  are  very  inefficient  at  doing  so  and  make  a  negligible 
contribution  [7]. 

Atmospheric  scatter  and  the  resulting  blur  introduced  to 
imaging  systems  has  been  well  researched  [8,  12],  but  there 
is  a  paucity  of  work  in  the  literature  looking  at  its  impact 
on  laser  dazzle.  The  authors  have  recently  developed  a  set  of 
simple  calculations  to  allow  the  visual  impact  of  laser  eye  dazzle 
to  be  modeled  [4] .  The  model  outputs  a  minimum  distance  for 
visual  detection  of  a  target  in  the  presence  of  laser  dazzle, 
known  as  the  nominal  ocular  dazzle  distance  (NODD). 
However,  the  initial  version  of  the  model  only  accounts  for 
intraocular  scatter  and  so  the  authors  are  seeking  to  incorporate 
the  contributions  of  extraocular  scatter  to  laser  eye  dazzle. 
Therefore,  this  experimental  paradigm  aims  to  understand 
the  relative  contribution  of  atmospheric  scatter  by  comparison 
to  the  CIE  general  disability  glare  equation  [13]  and  a  modified 
version  of  this  equation  [4] ,  both  of  which  estimate  the  scatter 
within  the  eye. 

2.  METHOD 
A.  Principles 

The  experiments  were  designed  to  derive  an  equation  for  the 
angular  distribution  of  atmospherically  scattered  laser  radiation, 
sometimes  referred  to  as  a  scatter  or  phase  function.  This  scatter 
function  describes  the  fraction  of  light  scattered,  per  unit  solid 
angle,  in  a  given  direction. 


Fig.  3.  Data  collection  schematic  showing  (top)  wide  acceptance 
angle  detector  (WAAD)  and  (bottom)  narrow  acceptance  angle  detec 
tor  (NAAD). 


Figure  3  shows  a  schematic  of  the  experimental  design. 
A  laser  was  propagated  over  a  distance  and  then  a  wide  accep 
tance  angle  detector  (WAAD),  located  on  axis  and  normal  to 
the  laser  path,  was  used  to  collect  light  from  all  angles  and  all 
directions  in  a  hemisphere  around  the  detector  plane.  The 
WAAD  gathered  both  on  axis  and  off  axis  forward  scattered 
radiation  from  the  laser  source  to  provide  a  measure  of  the  total 
irradiance,  I total (thW  •  cm"2).  A  narrow  acceptance  angle 
detector  (NAAD)  was  then  used  to  sample  around  this  hemi 
sphere  and  measure  the  contribution,  7(0) (mW  •  cm-2),  to 
I total  within  a  fixed  solid  angle,  Q  (sr),  at  a  given  angle,  6 
(deg),  from  the  laser  axis  in  a  horizontal  plane.  The  NAAD 
remained  within  the  center  of  the  laser  beam  but  was  rotated 
at  varying  angles,  replicating  the  way  an  eye  might  rotate  to 
view  a  scene  when  subjected  to  laser  dazzle. 

By  gathering  data  at  a  range  of  angles,  the  atmospheric  scat 
ter  function,  f{0)  (sr-1),  can  be  expressed  as  the  fraction  of 
light  scattered  per  unit  solid  angle: 

=  W/Q'  (1) 

f  TOTAL 

For  a  given  laser  irradiance  at  the  eye,  U (W  •  m'2),  the  re 
suiting  veiling  luminance,  Z*,(cd  ‘  m  2)>  can  then  be  calculated 
by  the  following  relationship  [14]: 


Lv  =  f{0)  •  683  -Vx-U,  (2) 

where  Vx  is  the  eye’s  photopic  efficiency  at  the  laser  wave 
length,  A  (0.883  for  332  nm),  and  683  is  the  maximum 
photopic  luminous  efficacy  in  lumens  per  watt  [15]. 


Research  Article 


Vol.  54,  No.  25  /  September  1  2015  /  Applied  Optics  7569 


B.  Implementation 

The  experiment  was  conducted  at  the  Outdoor  Directed 
Energy  Site  of  the  Tri  Service  Research  Laboratory  on  Joint 
Base  San  Antonio  Fort  Sam  Houston,  Texas.  This  outdoor 
laboratory  consists  of  a  400  m  gravel  track  with  a  small  building 
at  one  end  (the  “send”  point)  housing  the  laser,  aligned  with  a 
4  m  high  brick  backstop  at  the  other  end  (the  “receive”  point) 
where  the  measurement  equipment  was  located  inside  an  open 
trailer.  The  propagation  distance  from  laser  aperture  to  detector 
was  measured  as  380  m  using  both  a  rangefinder  and  a  global 
positioning  system  receiver. 

The  laser  source  was  a  commercially  available  332  nm  laser 
system  (Coherent  Verdi  V  18),  chosen  as  its  wavelength 
matched  the  most  commonly  used  laser  dazzle  devices,  while 
its  higher  power  (~15  W)  improved  the  ability  to  measure 
small  amounts  of  scatter.  The  optical  delivery  system  was  de 
signed  to  minimize  system  scatter  by  using  a  single  broadband 
alignment  mirror  to  direct  the  raw  beam  from  the  laser  aperture 
toward  the  receive  point  (Fig.  4).  During  nonmeasurement 
procedures  a  second  mirror,  mounted  to  a  motorized  flipper 
mount,  was  inserted  into  the  path  to  direct  the  laser  output 
to  a  beam  dump.  With  a  divergence  of  around  0.3  mrad, 
center  beam  irradiances  at  the  receive  point  were  approximately 
200  mW  •  cm“2. 

The  key  measurement  equipment  located  at  the  receive 
point  was  a  detector  system  comprising  two  laser  detectors 
(Thorlabs  PDA100A)  mounted  on  a  two  axis  translation  stage 
providing  vertical  and  horizontal  positioning  (Fig.  5).  The 
upper  detector  was  the  WAAD  with  a  ±50°  acceptance  angle 
to  give  the  /total  reading,  under  the  assumption  that  the  sig 
nal  beyond  this  acceptance  angle  was  negligible.  The  lower  de 
tector  was  the  NAAD,  with  a  designed  1°  acceptance  angle  and 
actual  value  of  0.955°  giving  a  solid  angle,  £2,  of  2.18  x  10"4  sr. 
The  NAAD  was  mounted  on  a  graduated  rotation  stage  such 
that  it  could  be  positioned  in  azimuth  independent  of  the 
upper  detector. 

The  NAAD  is  pictured  in  Fig.  6  and  consists  of  a  25.4  mm 
diameter,  60  mm  focal  length  achromatic  lens  to  focus 
incoming  light  through  a  precision  1  mm  aperture  in  front 
of  the  detector.  The  effective  input  aperture  diameter  was 
2.291  cm,  which  was  the  inner  diameter  of  the  retaining  ring 
used  to  hold  the  lens  in  place.  A  summary  of  the  NAAD  di 
mensions  is  given  in  Table  1.  The  detector  was  mounted  to  the 
rotation  stage  with  the  lens  positioned  at  the  center  of  rotation. 


Fig.  4.  Laser  optical  table  setup. 


Fig.  5.  WAAD  (top  aperture)  and  NAAD  (lower  aperture)  mounted 
on  translation  and  rotation  stages. 


Fig.  6.  Computer  aided  design  drawing  of  the  NAAD  showing  the 
lens  (right),  aperture  (middle),  and  detector  surface  (left).  The  inset 
shows  a  drawing  of  the  fully  assembled  NAAD. 


Table  1 .  Dimensions  and  Calculations  for  the  NAAD 


Parameter  Value 


Lens  focal  length  (mm) 

60 

Detector  aperture  diameter  (mm) 

1 

Calculated  field  of  view  (deg) 

0.955 

Calculated  field  of  view  (sr) 

2.18  x  10  4 

Input  aperture  diameter  (cm) 

2.291 

Input  aperture  area  (cm2) 

4.122 
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As  the  assembly  was  rotated  with  respect  to  the  direct  incoming 
beam,  the  light  focused  by  the  lens  translates  off  the  aperture 
blocking  it  from  reaching  the  detector.  Only  rays  within  the  1° 
acceptance  angle  reach  the  detector. 

Both  detectors  were  fitted  with  532  nm  laser  line  filters 
(Thorlabs  FL532  3)  to  enhance  their  sensitivity  by  isolating 
the  narrowband  laser  signal  from  the  broadband  ambient 
lighting.  In  their  fully  configured  states,  both  detectors  were 
calibrated  for  power  measurements  at  each  of  their  seven  gain 
settings  against  a  calibrated  detector  and  optical  power  meter 
(Newport  918D  and  2936  C).  For  the  NAAD,  the  calculated 
total  received  power  (W)  was  converted  to  irradiance 
(W  •  cm-2)  through  division  by  the  input  aperture  area  as  given 
in  Table  1.  For  the  WAAD  the  power  was  divided  by  the 
detector  active  area  (0.754  cm2)  to  give  irradiance. 

During  each  data  point  collection  both  detectors  operated  at 
their  native  1 0  kHz  rate  over  a  period  of  60  s  to  average  out  the 
effects  of  atmospheric  turbulence.  A  Lab  VIEW  software  inter 
face  allowed  the  operators  to  obtain  mean  and  standard 
deviation  figures  across  the  data  collection  period.  The  standard 
deviation  as  a  percentage  of  the  mean  /total  >  (StdDev%),  was 
used  as  an  approximate  measure  of  the  turbulence,  under  the 
assumption  that  other  atmospheric  processes  did  not  vary  dur 
ing  the  data  collection  period. 

Critical  to  the  scatter  measurement  system  was  the  inclusion 
of  an  obscuration  disc  in  front  of  the  NAAD,  designed  to  block 
direct  incident  rays  (see  Fig.  7  and  Fig.  8).  While  the  direct 
beam  was  not  within  the  field  of  view  of  the  NAAD,  the  ob 
scuration  disc  was  necessary  to  prevent  the  high  power  direct 
beam  from  entering  the  lens  of  the  NAAD  from  where  internal 
reflections  could  potentially  reach  the  detector  surface.  The  size 
of  the  disc  was  slightly  larger  than  the  NAAD  aperture  at 
approximately  25  mm,  which  facilitated  alignment  by  casting 
a  shadow  across  the  NAAD  aperture.  The  obscuration  disc  had 
to  be  located  greater  than  2.51  m  away  from  the  NAAD  aper 
ture  to  ensure  that  it  did  not  encroach  upon  that  detector’s  field 
of  view  and  block  any  scattered  light  that  was  part  of  the  desired 
measurement. 

Also  used  in  the  setup  was  a  camera  (Cohu  4815)  to  profile 
the  beam  prior  to  each  experiment.  The  camera  also  aided  in 
aligning  the  detectors  to  the  center  of  the  beam  via  the  trans 
lation  stages,  as  can  be  seen  in  Fig.  9. 

To  record  atmospheric  conditions  during  the  experiments,  a 
wireless  weather  station  (Davis  Vantage  Pro2  Plus)  was  located 
at  the  test  site.  Captured  data  of  relevance  included  tempera 
ture,  humidity,  atmospheric  pressure,  and  solar  radiation.  For 
atmospheric  quality,  the  concentration  of  particles  less  than 


Obscuration  disc  NAAD 


Fig.  7.  Illustration  of  the  obscuration  disc  used  to  block  the  direct 
beam  from  entering  the  NAAD  optics. 


Fig.  8.  Obscuration  disc  on  the  left  mounted  on  a  thin  pole,  casting 
a  shadow  across  the  aperture  of  the  NAAD  on  the  right  (close  up  lower 
right). 


Fig.  9.  Beam  profile  image  used  for  detector  alignment. 


2.5  pm  diameter  (PM  2.5  in  units  of  pg  •  m"3)  was  recorded 
based  upon  data  published  by  the  Texas  Commission  on 
Environmental  Quality  from  a  continuous  air  monitoring  sta 
tion  (CAMS)  located  approximately  four  miles  south  of  the  test 
site  (City  Public  Service  Energy  CAMS  678)  [16].  As  a  secon 
dary  measure  of  air  quality  the  meteorological  visibility  was  re 
corded  from  weather  data  at  San  Antonio  International  Airport 
[17],  which  is  located  approximately  five  miles  north  west  of 
the  test  site. 

C.  Procedure 

The  laser  was  first  directed  to  be  roughly  incident  upon  the 
detectors  by  adjustment  of  the  laser  output  alignment  mirror 
at  the  send  point,  facilitated  by  radio  communications  with  the 
receive  point.  The  detector  system’s  translation  stage  was  then 
moved  to  bring  the  WAAD  detector  into  the  center  of  the 
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beam  while  viewing  the  image  from  the  beam  profile  camera 
and  maximizing  the  received  laser  irradiance  reported  from  the 
detector. 

To  ensure  that  the  WAAD  detector  surface  was  normal  to 
the  laser  beam,  a  small  aperture  was  temporarily  placed  approx 
imately  1  m  in  front  of  the  detector  and  the  detector  was  ad 
justed  so  the  backreflection  from  the  detector  surface  showed  as 
a  bright  spot  on  the  rear  of  the  aperture.  The  alignment  of  the 
detector  was  then  fine  tuned  to  ensure  that  the  reflected  beam 
traveled  perfectly  back  through  the  aperture,  thus  indicating 
normal  alignment  with  the  laser  beam. 

Before  each  measurement  commenced,  the  detector  was 
zeroed  by  blocking  the  beam  at  the  source  and  initiating  a  zero 
ing  routine  in  the  software  to  remove  ambient  light  from  the 
readings.  After  unblocking  the  laser,  a  60  s  data  collection 
period  was  initiated  to  derive  the  initial  value  for  /total 
(typically  around  200  mW  •  cm'2,  as  previously  noted). 

The  NAAD  was  then  centered  in  the  beam  with  the  trans 
lation  stage  and  set  to  be  normal  to  the  laser  beam  using  the 
same  alignment  method  as  used  previously.  It  was  then  rotated 
about  its  center  axis  to  a  starting  angle  of  1°  and  the  obscuration 
disc  was  placed  at  the  predetermined  distance  of  greater  than 
2.51  m  from  the  detector.  The  obscuration  disc  was  first 
aligned  by  eye  to  ensure  that  its  shadow  completely  covered 
the  NAAD’s  aperture,  and  then  alignment  was  fine  tuned  by 
minimizing  the  measured  irradiance  on  the  detector.  The  start 
ing  angle  of  1  °  ensured  that  the  obscuration  disc  was  not  within 
the  detector’s  field  of  view,  and  also  that  the  Arago  spot  (the 
bright  spot  at  the  center  of  the  disc’s  shadow  caused  by 
Fresnel  diffraction)  had  no  measurable  impact  on  the  detector. 

At  this  point  the  NAAD  irradiance  was  zeroed  and  the 
60  s  data  collection  began  to  capture  the  value  for  /(1°),  which 
was  typically  around  200  nW  •  cm-2.  The  detector  was  sub 
sequently  rotated  by  0.5°,  following  which  the  obscuration  disc 
alignment  was  fine  tuned,  the  NAAD  was  zeroed  once  more 
and  data  gathering  began  again.  This  procedure  was  repeated 
until  the  signal  detected  by  the  NAAD  reached  its  sensitivity 
floor  of  approximately  1  nW  •  cm'2,  which  generally  occurred 
around  10°  from  the  laser  axis.  To  verify  that  the  small  signals  at 
large  angles  were  from  the  laser,  and  not  noise,  a  chopper 
wheel  was  used  at  the  laser  output  aperture  to  encode  a 
100  Hz  signal  on  the  beam.  Viewing  the  received  irradiance 
in  the  frequency  domain  showed  a  clear  peak  at  the  chopping 
frequency,  thus  verifying  that  the  signal  was  indeed  from  the 
laser  source. 


Finally,  the  NAAD  was  moved  back  to  1°  and  another  mea 
surement  taken  in  order  to  ensure  consistency  across  the  data 
collection  period,  which  typically  took  around  1  h.  The  WAAD 
detector  was  then  moved  back  into  the  center  of  the  beam  and 
another  value  for  /total  measured,  which  was  subsequently 
averaged  with  the  initial  measurement  to  give  the  final  value. 
As  can  be  seen  from  Fig.  9  the  beam  size  was  not  large  enough 
to  permit  continuous  monitoring  of  /total  throughout  the 
experiment,  and  so  the  averaging  of  start  and  end  readings  pro 
vides  some  compensation  for  any  changes  in  beam  quality,  tur 
bulence,  or  atmospheric  transmission  across  the  test  period. 

With  /total  and  aU  /($)  data  collected,  the  scatter  values, 
f(0),  could  then  be  calculated  using  Eq.  (1). 

3.  RESULTS  AND  DISCUSSION 

Data  were  collected  on  nine  separate  sessions  across  June  and 
July  2014.  A  summary  of  the  atmospheric  conditions  for  these 
dates  is  given  in  Table  2.  Atmospheric  parameters  were  aver 
aged  over  the  data  collection  periods  (approximately  1  h  each) 
wherever  possible,  or  taken  from  time  periods  in  the  middle  of 
the  collection  period. 

Testing  spanned  dates  with  good  (PM  2.5  less  than 
12  pg  •  m“3)  and  moderate  (between  12  and  35  pg  •  m“3) 
air  quality  [16],  with  temperatures  ranging  from  68°F  to 
95°F  and  humidity  ranging  from  51%  to  93%.  Two  data  ses 
sions  took  place  at  night  with  no  measurable  solar  radiation, 
while  the  others  took  place  during  the  morning  or  peak  after 
noon  sunshine.  The  StdDev%  factor,  giving  a  rough  indication 
of  turbulence  from  the  standard  deviation  of  the  mean  /total  > 
ranged  from  12%  (low  turbulence)  up  to  75%  (high  turbu 
lence).  The  visibility  data  from  the  airport  location  usually  fol 
lowed  the  trend  of  the  PM  2.5  data,  with  concentrations  of 
more  than  14  pg  •  m“3  corresponding  with  visibility  less  than 
10  miles.  The  one  exception  was  data  for  09  July,  which  indi 
cated  the  highest  recorded  PM  2.5  concentration  but  a  visibil 
ity  of  greater  than  10  miles.  This  could  be  due  to  the  separation 
of  the  airport  from  the  CAMS  by  approximately  8.5  miles,  with 
the  possibility  of  localized  weather  conditions. 

An  aggregated  scatter  graph  across  all  nine  data  collection 
periods  is  shown  in  Fig.  10.  It  can  be  seen  that  the  atmospheric 
scatter  data  were  very  similar  for  each  collection  period,  with 
the  =t  1  standard  deviation  error  bars  (not  shown  on  this  graph 
for  clarity)  for  each  collection  period  mostly  encompassing  each 
other.  There  are  some  outliers  in  the  dataset,  most  notably  the 


Table  2.  Summary  of  Atmospheric  Conditions  for  the  Nine  Days  of  Data  Collection 


Date 

Time 

Solar 
(W  •  nr2) 

Temp. 

(°F) 

Humidity 

(%) 

Pressure 

(mb) 

PM-2.5 
(fig  ■  m"3) 

Visibility 

(miles) 

StdDev% 

(%) 

04  June 

09:00 

338 

78 

81 

1011.3 

6.3 

>10 

38 

05  June 

09:00 

280 

79 

79 

1009.1 

9.4 

>10 

30 

06  June 

09:00 

374 

80 

77 

1011.2 

25.3 

9 

60 

09  June 

09:00 

133 

78 

81 

1007.1 

14.4 

8 

21 

10  June 

05:00 

<1 

68 

93 

1008.9 

7.5 

>10 

12 

12  June 

15:00 

780 

95 

51 

1007.2 

13.4 

>10 

75 

24  June 

10:00 

454 

83 

75 

1015.0 

11.1 

>10 

61 

07  July 

09:30 

537 

81 

67 

1014.9 

25.4 

8 

69 

09  July 

04:30 

<1 

75 

86 

1014.7 

25.9 

>10 

17 
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Fig.  10.  Scatter  versus  angle  from  all  nine  data  collection  sessions. 


lowest  set  of  data  points  at  2.5°  to  4.5°  that  were  all  from  04 
June.  While  the  PM  2.5  for  that  day  was  the  lowest  of  the  nine 
sessions,  it  was  only  marginally  better  than  other  days  and  there 
were  no  other  noticeable  differences  in  the  atmospheric  oondi 
tions  for  that  day.  The  dissimilarity  with  the  subsequent  eight 
data  sessions  leads  the  authors  to  believe  that  an  experimental 
error  may  have  occurred  during  the  data  collection  possibly 
an  unintentional  movement  of  the  detector  bqrond  the  angle 
displayed  on  its  rotation  stage.  Data  from  04  June  were  there 
fore  removed  from  the  dataset  for  the  subsequent  analysis. 

Fig.  11  shows  a  plot  of  the  highest  (09  July)  and  lowest  (24 
June)  recorded  scatter  Sanctions  from  the  reduced  set  of  eight 
data  collections  sessions,  with  error  bars  indicating  rirl  standard 
deviation.  The  short  error  bars  for  09  July  indicate  low  turbu 
lence  (StdDev%  =  17%)  compared  with  the  large  error  bars  for 
24  June  (StdDev%  =  61%).  Within  the  measured  errors,  both 
scatter  functions  are  very  dose.  Looking  at  the  summary 
weather  statistics  it  could  be  argued  that  any  scatter  differences 
are  due  to  the  high  PM  2.5  concentration  on  09  July 
(25.9  pg  ■  m"3)  compared  with  24  June  (11.1  pg- nT3), 
However,  Fig.  12  shows  a  comparison  of  the  scatter  functions 
recorded  for  another  set  of  high  (07  July*  25.4  pg  -  m~3)  and 
low  (10  June,  7-5  pg  ■  m_?)  PM  2.5  values  where  there  is 
no  noticeable  difference  in  the  scatter  functions. 


Angle  (deg) 

Fig.  11.  Comparison  of  highest  and  lowest  scatter  captured. 


Angle  (deg) 

Fig.  1 2.  Comparison  of  scatter  for  highest  and  lowest  PM  2.5  con 
centrations. 


Fig.  13  shows  a  comparison  of  scatter  data  from  1 0  June  and 
12  June:  the  days  of  lowest  and  highest  solar  radiation  (<1  ver 
sus  780  W  ■  m~2),  lowest  and  highest  temperature  (68°F  versus 
95°F),  highest  and  lowest  humidity  (93%  versus  51%)*  and 
lowest  and  highest  StdDev%  turbulence  (12%  versus  75%). 
The  higher  turbulence  accompanies  the  high  temperature, 
high  solar  radiation  condition  to  give  much  Iaiger  uncertainties 
for  the  12  June  data,  but  the  mean  values  look  broadly  similar 
to  10  June.  The  similarity  given  the  vastly  different  solar  radi 
ation  levels  (dark  night  versus  peak  summer  day)  shows  that  the 
experimental  implementation  and  procedure  were  robust  to  in 
fluence  from  ambient  light.  It  also  indicates  that  the  measures 
taken  to  average  out  turbulence  effects  were  successful,  in  terms 
of  scatter  dependencies,  the  data  suggest  that  temperature  and 
humidity  do  not  affect  atmospheric  scattering,  at  least  across 
the  range  of  values  that  these  two  datasets  represent  and  when 
averaged  over  60  seconds. 

Fig.  14  shows  all  mean  data  from  the  reduced  set  of  eight 
collection  sessions,  together  with  the  following  line  of  best  fit 
scatter  function  (R2  =  0.95): 

f{0)  =  0.0031  (3) 


Angle  (deg) 

Fl  g .  1 3.  Co  mpaifco  n  of  scatter  fo  r  highest  and  1  owest  sol  ar  radiation 
and  temperature. 
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Angle  (deg) 

Fig.  14.  All  data  from  reduced  set  of  eight  data  collection  sessions 
together  with  line  of  best  fit. 


in  Fig.  I  5,  this  scatter  function  is  compared  with  the  glare 
due  to  scatter  in  the  human  eye  in  order  to  understand  the 
overall  contribution  of  atmospheric  scatter  to  laser  eye  dazzle. 
The  upper  curve  is  generated  from  the  CIE  general  disability 
glare  equation  an  internationally  accepted  standard  equation 
for  intraocular  scatter,  based  upon  human  subject  measure 
ments  [13].  The  CIE  equation  accounts  for  increasing  scatter 
with  age  and  lighter  eyes,  and  so  the  plotted  graph  is  for  the 
lowest  possible  scatter  case  of  a  20  year  old  with  very  dark  eyes. 
It  can  be  seen  that  scatter  from  the  eye,  even  in  this  lowest 
scatter  case,  is  still  over  3  orders  of  magnitude  greater  than 
the  scatter  from  the  atmosphere  found  in  the  present 
experiment. 

The  authors  have  previously  presented  an  adjustment  to  the 
CIE  general  disability  glare  equation  which  suggests  a  lower 
scatter  in  low  ambient  luminance  levels  wrhen  simulating  laser 
exposures  [4,18].  Using  this  adjustment  factor  with  a  night 
time  ambient  luminance  level  of  0.001  cd  -  m'2  gives  the 
middle  dotted  line  in  Fig.  15.  This  now  represents  the  lowest 
possible  predicted  intraocular  scatter,  based  on  a  20  year  old 
with  dark  eyes  viewing  in  the  lowest  ambient  lighting 


Fig.  15.  Comparison  of  the  atmospheric  scatter  measured  in  this 
experiment  to  the  predicted  eye  scatter  from  die  CIE  general  disability 
glare  equation,  and  a  version  adjusted  for  low  ambient  light  level  laser 
exposures. 


conditions.  Even  in  this  case,  the  predicted  atmospheric  scatter 
is  less  than  5%  of  the  predicted  eye  scatter. 

Referring  back  to  Eq.  (2),  each  scatter  function  can  be  used 
to  estimate  an  equivalent  wiling  luminance  distribution  on  the 
retina.  Taking  a  532  nm  laser  exposure  of  500  pW  ■  cm'2  (half 
the  maximum  permissible  exposure  for  a  10  s  visible  wave 
length  exposure,  based  upon  laser  safety  standards  [19])*  and 
an  angle  from  the  laser  axis  of  5°,  gives  an  equivalent  veiling 
luminance  caused  by  the  atmospheric  scatter  of  0.3  cd  -  m'2. 
This  compares  with  7  cd  -  m  2  caused  by  the  adjusted  low  light 
eye  scatter  case,  and  852  cd  ■  m"2  by  the  CIE  case.  The  veiling 
luminance  equation,  Eq.  (2),  maintains  the  relatiw  magnitude 
of  the  atmospheric  contribution,  with  it  being  less  than  5%  of 
the  luminance  caused  by  the  adjusted  eye  equation  and  less 
than  0.05%  of  the  luminance  of  the  CIE  q^e  equation. 

4.  CONCLUSIONS 

The  results  generated  by  this  study  suggest  that  a  robust  experi 
ment  has  been  devised  for  the  measurement  of  atmospheric 
scatter.  It  has  been  demonstrated  that  very  small  levels  of  atmos 
pheric  scatter  (^nW  -  cm"2)  can  be  measured,  independent  of 
ambient  light  level  and  atmospheric  turbulence. 

For  short  range  laser  engagements  in  atmospheres  of  good  to 
moderate  air  quality,  this  study  suggests  that  atmospheric  scat 
ter  has  a  negligible  impact  upon  laser  eye  dazzle.  Even  with  very 
conservative  estimates  of  the  scatter  contribution  from  the 
human  eye,  the  atmospheric  scatter  adds  no  more  than  5% 
to  the  overall  scatter  for  a  laser  propagating  over  a  distance 
of  380  m.  Such  scatter  levels  would  not  create  any  significant 
additional  veiling  luminance  across  the  scene  for  a  human 
observer  when  compared  to  the  veiling  luminance  introduced 
hy  the  scatter  within  the  eye. 

The  magnitude  of  atmospheric  scatter  does  not  show  any 
clear  correlation  with  atmospheric  quality  based  on  PM  2.5 
concentrations,  at  least  for  the  good  and  moderate  PM  2;5 
levels  experienced  across  the  380  m  range.  Higher  PM  2.5  con 
cent  rations  could  cause  an  increase  in  the  significance  of 
atmospheric  scatter,  although  the  present  experiment  cannot 
further  inform  such  predictions. 

A  statistical  study  of  historical  PM  2.5  data  from  the  CAMS 
678  site  [16]  shows  that,  from  2011  to  2013,  there  has  only 
been  an  average  of  43  hours  per  yrear  where  the  PM  2.5  con 
centration  has  risen  above  35  pg  ■  m"3  and  into  the  ^unhealthy" 
category.  This  generally  occurs  in  very  narrow  time  periods,  on 
average  for  3  hours  out  of  14  days  per  year.  Therefore,  the  re 
suits  presented  herein  can  be  applied  99-5%  of  the  time  during 
an  average  year  at  the  test  site:  To  explore  PM  2.5  scatter 
dependencies  further  would  require  relocation  of  the  experi 
ment  to  an  environment  with  much  poorer  air  quality  than 
generally  experienced  in  San  Antonio.  Alternatively,  the  use 
of  artificial,  localized  atmospheric  disturbances  (e.g.,  smoke) 
could  be  considered. 

It  is  anticipated  that  longer  propagation  ranges  would  in 
crease  the  significance  of  atmospheric  scatter,  due  to  the  greater 
volume  of  atmosphere  for  light  interaction.  Again,  the  present 
experiment  cannot  directly  inform  such  predictions,  but  the 
relocation  of  the  experiment  to  a  longer  range  wTould  provide 
additional  data  for  comparison.  The  use  of  large  mirrors  at  the 
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receive  point  of  the  current  test  site,  to  reflect  the  beam  and 
double  the  path  length,  is  another  option. 

Within  the  environmental  parameters  encompassed  by  this 
experiment,  atmospheric  scatter  appeared  to  be  independent  of 
temperature,  humidity,  and  atmospheric  pressure  given  equiv 
alent  PM  2.5  concentrations.  However,  the  collection  of  addi 
tional  data  during,  for  example,  very  cold  conditions  would 
provide  a  useful  comparison  to  further  confirm  this  finding. 

Additional  improvements  are  possible  to  better  understand 
how  air  quality  influences  this  measurement.  Having  the 
capability  to  measure  air  quality  at  the  test  site  would  be  pref 
erable  to  account  for  any  localized  variations,  and  this  would 
also  allow  artificial  pollutants  to  be  used  and  monitored.  Use  of 
a  scintillometer  to  measure  the  refractive  index  structure 
parameter,  Cn 2,  would  be  beneficial  to  characterize  turbulence. 
Using  a  higher  power  laser  would  also  provide  additional  signal 
for  measurements,  potentially  allowing  data  collection  beyond 
the  10°  limit  found  with  the  15  W  laser. 

Finally,  a  natural  extension  to  this  experimental  work  would 
be  the  introduction  of  a  theoretical  component.  Previously  re 
ported  models  [8  1 2]  could  be  validated  with  the  experimental 
data,  and  attempts  made  to  apply  the  simulations  to  longer 
ranges  and  higher  PM  2.5  concentrations. 
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